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Abstract. The complex interactions of localized vortices with waves is investigated using a model of point 
vortices in the presence of a transverse or longitudinal wave. This simple model shows a rich dynamical 
behavior including oscillations of a dipole, splitting and merging of two like-circulation vortices, and chaos. 
The analytical and numerical results of this model have been found to predict under certain conditions, the 
behavior of more complex systems, such as the vortices of the Charney-Hasegawa-Mima equation, where 
the presence of waves strongly affects the evolution of large coherent structures. 
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Physical systems like large scale motion in oceans and at- 
mospheres [1,2 ,3J, space and laboratory plasmas [4,5j, or 
in statistical physics, the XY planar spin model [6,7J, are 
all able to generate both vortices and waves, which interac- 
tion is one of the keys to their behavior. The phenomenol- 
ogy of these systems is very rich, ranging from turbulence 
to self-organization, formation of coherent structures ^,9^ 
[Tm fTT] and anomalous transport [T2l[T3l[T3] . Understand- 
ing the elementary mechanisms dominating the dynamics 
is crucial to model the statistical properties of these flows. 
One important fundamental process is the vortex interac- 
tion in presence of waves. Besides complex dynamical pat- 
terns, particular behaviors have been identified, indeed, it 
has been observed that collisions of vortices may induce 
their merging, and on the other hand, the inverse process 
of splitting is also possible, large vortices may break under 
the action of waves. Both mechanisms are accompanied by 
the emission of waves, as for instance in the interaction of 
plasma vortices and drift waves [T5j[T6l[T7] . 

The role played by waves in the interaction and evo- 
lution of vortices is not well understood. A comparison of 
the long-time evolution of Euler decaying two-dimensional 



servation equation which is a slight generalization of the 
vorticity Euler equation. In general the existence of vor- 
tex solutions is related to the so-called Poisson nonlinear- 
ity, which is the two-dimensional version of the convective 
nonlinearity in a normal fluid, or the electrostatic drift 
in magnetized plasmas. Besides, waves may be related to 
some source of divergence of the velocity field, as com- 
pressibility or as in the two-dimensional reduction of a 
stratified flow. The equation for the generalized vorticity 
i7, containing these basic mechanisms (a linear dispersion 
relation, and a Poisson nonlinearity) is given by 



Si? 



(1) 



where [/, g\ = f^gy - fyg^ (/ and g are arbitrary func- 
tions of and the subscripts denote derivatives) is 
the usual Poisson bracket, and tjj is the stream function. 
The actual relation f2 = F{tlj) depends on the considered 
physical system, for the Euler equation it is simply given 
by f2 = — V^V^. This equation expresses the conservation 
of the generalized vorticity following the current lines. The 
Hasegawa-Mima equation for plasma drift waves [41122]. 
also known in fluid mechanics as the Charney geostrophic 



turbulence [lOj and the Charney-Hasegawa-Mima one (which equation for Rossby waves i23,.24,25i, reduces to with 
supports wave modes) [H], shows that the relaxation to- the generalized vorticity. 



wards a kind of thermodynamic equilibrium state [1^] is 
much longer or unreachable in presence of a wave field [T^ 
[2Qt [2T] . This illustrates that vortex interactions are largely 
influenced by the presence of waves. 

A common feature of these vortex- wave systems is that 
they are two-dimensional and may be described by a con- 



(2) 
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where is, in the plasma case, related to the electric po- 
tential(in suitable units), ps is the hybrid Larmor radius 
(Rossby length in the atmosphere) and Vd is the drift wave 
velocity. The similarity of the equations governing Rossby 
waves in the atmosphere and drift waves in a plasma. 
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may be traced back to the analogy between Coriolis and 
Lorentz forces, which determine the form of the dispersion 
relation, and the common convective nonlinearity, which 
reduces to a Poisson bracket in two dimensions. These two 
parameters, a length scale and a characteristic velocity, are 
related to important physical effects that modify the sys- 
tem's dynamics with respect to the one of the simpler Eu- 
ler fluid. Indeed, the Charney-Hasegawa-Mima equation 
admits both waves and localized vortices solutions fT7l[26l 
[77|[^ [TT] , allowing a rich nonlinear dynamics as shown in 
numerical simulations. For instance, a dipole vortex has an 
oscillatory trajectory when its symmetry axis is inclined 
with respect to the wave direction of propagation [29], or 
like-sign vortices (monopoles) may merge, under suitable 
conditions JT5\. 

The drift velocity is related to the phase velocity of 
drift waves, the dispersion relation being 

= Vdky/(l + plk^) , (3) 

with uj the frequency and k the two-dimensional wavenum- 
ber. In the limit of vanishing drift velocity, localized vor- 
tices can become point vortices, 

Q ^ ^5{x) (4) 

of circulation i^, with an interaction range of the order of 
the characteristic length pg. 

The vortex-wave interaction is in general extremely 
complicated, involving the emission and scattering of waves 
by vortices and simultaneously the deformation of the vor- 
tex shape by the wave field [30]. In this paper we investi- 
gate the interaction of localized vortices with waves, using 
the approximation that the wave is fixed, not affected by 
the vortices. We distinguish two situations according to 
the longitudinal or transverse character of the imposed 
wave. Transverse waves are naturally generated by the 
dynamics of Eq. ([T]), while longitudinal ones appear as 
a source term associated with the longitudinal component 
of the pressure driven velocity [31,32j. We show that in 
both cases the interaction of vortices (with same topologi- 
cal charge) with the wave can trigger their merging, or de- 
pending on the initial positions and wave strength, their 
separation. Moreover complex, chaotic or quasiper iodic, 
behavior can be found. We note as well that event houg in 
a different context, similar features can be observed when 
"perturbing" the integrable motion of two equal-sign vor- 
tices with noise [33j[3l]. 

In the following section, Sec.[2l we state the basic equa- 
tions of the point vortex model for the transverse and 
longitudinal external wave, and discuss its relation with 
Charney-Hasegawa-Mima equation. In Sees. [3] and H] we 
investigate the different dynamical regimes that the pres- 
ence of the wave, transverse and longitudinal respectively, 
can induce in the motion of the point vortices. We de- 
scribe the dipole (opposite circulation vortices) and the 
monopole (like-circulation vortices) trajectories. We find 
the range of parameters where the strong interaction of the 
monopole and the wave leads to their merging or splitting. 



In Sec.[5]we compare the results of Sees. [3] and HI with nu- 
merical simulations of the Charney-Hasegawa-Mima equa- 
tion, and finally (Sec. [6]) we conclude with a summary and 
brief discussion of the results. 



2 Transverse and longitudinal waves: point 
vortex model 

In principle, the wave term in the generalized vorticity ([2]), 
being proportional to the x-coordinate and obviously non 
localized, forbids a solution in the form of an assembly of 
point vortices [35j. However, as a first approximation, we 
may consider a system where the vorticity is highly con- 
centrated but subject to the action of a wave field in such 
a way that the effect of the vortices on the wave may be 
neglected. If one neglects the internal structure of the lo- 
calized vortices and concomitantly, the feedback effects on 
the wave, the velocity field can be considered as resulting 
from the superposition of two terms, one related to the 
point vortices and one related to the wave. Therefore, our 
basic model will be two point vortices driven by the ac- 
tion of an external wave. It mimics the phenomenology of 
the Charney-Hasegawa-Mima system, dominated by the 
dynamics of coherent structures evolving in the field of 
transverse and driven longitudinal waves. We show that, 
in spite of its extreme simplicity, this system contains a 
rich dynamics that can model processes such as splitting 
and fusion, and deterministic or chaotic propagation of 
vortices. 

Point vortices are defined by a potential vorticity dis- 
tribution given by a superposition of Dirac functions, 

1 ^ 

^i{x,t) = — Y.r^S{x-Xo^{t)) , (5) 

where is a vector in the plane of the flow, is the 
circulation of vortex a, TV is the total number of vortices, 
and Xa{t) is the vortex position at time t. Point vortices 
are known to be exact solutions of Euler equation [36j, 
that is when the vorticity is given by 

f2 = -V^i'i , (6) 

in terms of the local stream function tpi , but they are also 
exact solutions of the more general equation ([T]), for 

f2 = f2i = -\/^^i^^i/pl. (7) 

The usual tjji (solution of the Poisson equation), which 
gives a logarithmic interaction, is replaced in the second 
case ([7j), by a modified Bessel function Kq interaction (so- 
lution of the Helmholtz equation). 

The presence of the wave should add to (j5]) a regular 
source of potential vorticity. Considering that the corre- 
sponding stream function is the superposition of a local- 
ized (vortex) and a wave term, = ?/;^ + ?/;^ , the motion 
of vortices will be governed by, 

Xa = v{x,t)\x^Xo.^ V = Z X V(V^^ + I/Jw) + V0 (8) 
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where the velocity field is calculated at the position of 
the vortex a, excluding the singular contribution of the 
a vortex and z denotes the unit vector parallel to the 

axis. In addition to the transverse wave, the last term 
0, takes into account the possibility of a longitudinal wave 
driven by an external source. 

In order to explicitly display the approximation un- 
derlying the point vortex model of Eq. (|8]), we return to 
the Charney-Hasegawa-Mima equation ([T]), that can be 
written as 

where f2i given by Eq. ([7j), is the localized potential vor- 
ticity. Introducing the decomposition tj; = tpy -\- ^1)^^ one 
obtains, 

d 

at Pi dy 

+^^-^- + [^-'-^'^-] = 0- 

pI dy 

(10) 

We consider now that the localized component of the po- 
tential vorticity, Qi = f2i{ipy)^ is given by the superpo- 
sition of point vortices ([5]). Then, the first line of ^TO\i 
only contains ^-localized terms, and vanishes if the point 
vortices trajectories are precisely given by (|8|) (without 
the source term, in V0). Therefore, the point vortex- wave 
model requires the fulfilling of the following two condi- 
tions. First, it assumes that the terms on the last line 
of ([To]), that is the dispersion effect on the vortices and 
the influence of the vortex on the wave, are negligible. 
Second, we remark that for a single drift wave (a wave 
satisfying the linear dispersion relation) the second line 
vanishes identically. 

In general ([9]) is complemented with a source and dis- 
sipation terms, that we write as, 

+ + [V^, = z^vV + s{x, t) , (11) 

dt pI dy 

where v is the kinematic viscosity, and S(x^t) an exter- 
nal vorticity source. In the Appendix [3 we show that the 
source term can be naturally associated with longitudinal 
wave perturbations. 

In summary, we apply the velocity field derived from 
the stream function whose dynamics is governed by (|9]) or 
(pTj) . to the point vortex model ([8]). In the source free case 
([9]), the velocity field is purely transverse, and we may 
write it as, 

v{x,t) = z xVil^^x.t) , t/j = i/ji ^ t/jyj (12) 

where the stream function is assumed to result from the 
superposition of the point vortex term ipi and the non- 
local (wave) term t/j^j. In contrast, in the forced case the 
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Fig. 1. (Color online) Poincare map of the two vortices- wave 
system (p^ , in {x = xi — X2,y = yi — y2) phase plane (a); run- 
away trajectories (vortex splitting) are present together with 
others approaching the origin (vortex merging) , and large chaos 
regions, (b) Zoom of the central region that show "merging" tra- 
jectories, an ensemble of initial points (blue circle) approaches 
the center (red arc) in one period. 

velocity field may not be transverse, and contain a longi- 
tudinal (potential) contribution. For the driven case ([TT]) . 
the velocity field can then be chosen as, 

v{x, t) = zx V^Iji{x, t) + V^(cc, t) , (13) 

in accordance with (|^ but with ^-independent quanti- 
ties, where, in addition to the rotational component de- 
rived from the stream function ?/^|, there is a gradient 
component, whose physical origin is related to the non- 
stationary pressure perturbations in the vertical direction 
[24[[3T] [32) . The first case ([12]) corresponds to a system of 
point vortices in the field of a transverse wave, and the 
second case ([T3|) . to the presence of a longitudinal wave. 

The motion of the vortices, as a consequence of the 
Helmholtz theorem [37J, is determined by the value of the 
velocity v = z x S/t/j + [which contains both cases 
([T2|) and ([T3|) ]. at the position of the vortices, as already 
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shown in (jSj). If the driven potential (j) vanishes, this mo- 
tion is Hamiltonian. In the general case the point vortex 
equations can be written as. 



dH 

dXa 



r,V(/)(x«,t). (14) 



where the Hamiltonian H is given by 

H = ^Y1 ^c.r^u{\xa - x^\) + , (15) 



27r , 

a>/3 a 

with U{x) = — log(x) in the Euler case, and U (x) = 
Ko(x/ps), the modified Bessel function of zero order, in 
the more general case (when ^ oo it tends to the loga- 
rithm). When the distance between the vortices is smaller 
than the typical interaction length ps the behavior of the 
two systems is similar, in the opposite case, the Kq in- 
teraction decreases exponentially and the vortices become 
"almost free". 

The Hamiltonian (p!5|) explicitly depends on time, and 
in the case of two vortices, reduces to the class of one and 
a half degrees of freedom, which generically display Hamil- 
tonian chaos [38j. We can therefore expect that two ini- 
tially well separated trajectories may approach each other 
closely or, at variance, two initially close vortices may fol- 
low diverging trajectories. Translating these behaviors to 
the extended system described by the Charney-Hasegawa- 
Mima equation, we anticipate merging and splitting of lo- 
calized structures. 



3 Dynamics of two point vortices in the field 
of a transverse wave 

We assume that a particular solution of (O, as discussed 
in the previous section, can be approximated hy i/j = i/ji -\- 
V^^y, where the first term represents two point vortices of 
circulations Fi and / 2, and the second term, a wave t/j^ = 
Fq cos{ky — uot) of amplitude /q, wavenumber k = (0, /c), 
and frequency uo. The motion equation ([T^ of the vortex 
pair writes. 



X2 



27rps V^2 -xi J r 
A fy2-yi\^i{r/Ps) 



2ttPs \Xl-X2 



(16a) 

xF^k sm{ky2 — oot) , 
(16b) 



the relative motion [31]. Using the coordinates {x^y) = 
{xi -X2,yi-y2)/2 and (X, F) = (xi +X2,^i +^2)/2, the 
relative motion is described by the set, 

X = — I — I - Fok cos( kY -ut) sin ky, (17a) 
27r psV \psj 



y 



(17b) 



where r = 2(x^ + y'^Y^'^ and Y = const. The important 
point about these equations is that they can be derived 
from an effective Hamiltonian: 

Ht = — Ko f — 1 Focos( kY -ujt) cos ky (18) 
27r \psj 

which describes the one and a half degrees of freedom 
dynamics in the (x, y) phase space, and constitute then 
a convenient representation to introduce a Poincare map, 
using the natural period 27r/cj of the driven wave. 

To identify the diverse dynamical patterns of ([T6|) , cor- 
responding to a set of initial conditions, we represent in 
Fig. [1] the Poincare map of the phase space r = (x, y) 
using the parameters F = Stt, Fq = 5.6, k = Stt/IO and 
uj = k/{l-\-k'^) (in units of ps and Vd). We may distinguish 
periodic and quasiperiodic, chaotic and runaway trajec- 
tories, depending on the initial positions of the vortices. 
Chaotic trajectories may include a set of initial conditions 
that rapidly (in a few periods) approaches the origin, as 
shown in Fig. [IJb). This particular kind of trajectories, 
in an extended dissipative system, should correspond to 
merging of vortices. Moreover, runaway trajectories should 
correspond to splitting of initially bounded monopoles. 
We show in Sec. [5l that these processes, including com- 
plex monopole interactions, are effectively observed in the 
dynamics of the Charney-Hasegawa-Mima equation. 

If the two vortices have opposite circulations Fi = 
— /2, the system ([T6]) becomes integrable. 



X = —Fok sin ky cos{kY — uot) , 
^ = 0, 

Fy Ki{r/ps) 



(19a) 
(19b) 



X = 
Y = 



2iTps r 
Fx Ki{r/ps) 
2iTps r 



Fok cos ky sm{kY — out) , (19c) 
(19d) 



where Xi = (xi^yi)^ z = 1, 2, are the vortex positions, r 
their distance of separation, Ki the modified Bessel func- 
tion of first order, and x the unit vector parallel to the 
X— axis. In the absence of wave, the system behaves as a 
nonlinear oscillator whose frequency is (A -{-F2 ) Ki (ro /ps ) /4 
(where tq is the initial distance between the two vortices). 
However, driven by the wave, the dynamics of this oscil- 
lator will generically exhibit chaotic trajectories. 

In the case of two identical point vortices of circulation 
F = Fi = F2, the center of mass motion separates from 



from which one finds that y = const., and r = r{t) sep- 
arates from the center of mass equations. Hence, the so- 
lutions of these equations are typically dipole trajectories 
that follow a straight line to which is superposed an os- 
cillation (Fig. [2^). One particular solution is obtained in 
the case y = n7r/k with n integer, the vortex separation 
in the y direction is a multiple of the wavelength. In this 
case the vortices drift with a constant velocity in the y 
^r^^tion and execute horizontal excursions in the form, 

X{t) = Vot - Fok cos{t/to) , (20) 

where the constants Vq and to depend on the initial condi- 
tions. These horizontal excursions are present in the tra- 
jectory shown in Fig. [2]3. 
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Fig. 2. (Color online) Trajectory of the dipole vortex in a 
transverse wave. Parameters: g = Stt, a = 4, xi — (0, 0) (a) 
A representative case of the oscillatory dipole motion X2 = 
(1, —0.5); (b) particular solution showing horizontal excursions 

X2 (0, -Tt). 



4 Dynamics of two point vortices in the field 
of a longitudinal wave 

We consider now the motion of two vortices in the presence 
of a longitudinal wave. The basic equations of the model 
are directly derived from ([HI), where we put ip^j = and 
choose a wave potential simply defined by = /q cos{k ■ 
X — out), with Fq the wave amplitude, k the wave vector, 
and the wave frequency. 



r = r^z X — — 2/o/c cos {k R — uot) sin f 



R = r_zx 



2r2 



- Fok sin (k • R — ujt) cos 



k • r 



, (21) 



, (22) 



(in the ps ^ oo limit), where r = ri — r2 is the relative 
position of the vortices, R = (ri + r2)/2 is their "cen- 
ter of mass", and r+ = (A + r2)/27r, r_ = (A - r2)/27r, 
are their reduced circulations. For the analytical computa- 
tions we take the simpler logarithmic interaction, although 
generalization to the modified Bessel interaction, should 
be straightforward. 

For further calculations, it is useful to measure lengths 
in 1/k and time in l/cu units, and to perform a Galilean 




X 



(b) 



Fig. 3. (Color online) Trajectory of the dipole vortex in a lon- 
gitudinal wave, (a) A representative case of the dipole motion 
locking when \V\ < ao; (b) the weak wave case (|y| > ao), the 
dipole trajectory is modulated by the wave. Parameters: g = I, 
a = 0.3, xi = (0, 0), X2 = (1, -0.7) (a) and X2 = (1, -1) (b), 
tf = 50. 



transformation to the frame moving at the wave phase ve- 
locity, leading to the following non-dimensional variables 
and parameters: = k - y' = k • {z x r)/2, f = ujt^ 
X' = k-R-t,Y' = k'{zxR),i\ie distance r' = (x', ^'), 
and center of mass R = (X',y) vectors, and parame- 
ters ao = /o/c^/cj, for the wave amplitude = F^k'^ /4uj, 
O- = F-k!^ /4(jj^ for the vortices circulation (in the follow- 
ing we drop the primes). The coordinate system is chosen, 
without loss of generality, so that the x-axis is parallel to 
the vector k. 

In addition to the analytical calculations, we compute 
the point vortices trajectories in the presence of a longitu- 
dinal wave, by direct integration of the motion equations 
in the form, 

i2 = ( 5 : l\ ) - sin(.. - t) , (23b) 

for arbitrary v = l/kpg. We used g = kF/uops^ where F = 
Fi = —F2 for the dipole (minus sign), and F = Fi = F2 
for the monopole (plus sign), and a = kFo/ujpg (in units 
of k = 1 and u = 1). The vortex motions are qualitatively 
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X— component of the speed of the unperturbed dipole in 
the frame of the wave. We obtain for R 



Fig. 4. (Color online) Quasi-periodic (a) and chaotic (b) 
motion of the monopole, in the long-wave case. Parameters 
xi = (0,0), X2 = (0.3,0) (a), and X2 = (1,0) (b), tf - 200, 
Ti = 1, = 1, To = 0.5, cj = 1. 



similar for finite range and logarithmic interactions; we 
show in the following, results for = 1, and different 
choices of the initial positions and strengths g and a. 

Two cases are of special interest, the so called "dipole", 
when the vortices have opposite circulations = 0, and 
the "monopole", when the vortices have the same circula- 
tion a_ = 0. Without wave the dipole moves in a straight 
line, and the monopole follows a circular trajectory. 

In the dipole case, the equations of motion reduce to 
the completely integrable system 



X = —ao cos X sin X , 
^ = 0, 

y 

X = — 1 — ao sin X cos x + o- — , 
Y = — a_— . 



(24a) 
(24b) 

(24c) 
(24d) 



where = ^ y'^ . Examples of dipole trajectories are 
presented in Fig. [3l 

Although the motion of a dipole interacting with a 
wave is integrable, to find a quadrature is difficult, so some 
insight can be gained by studying particular exact solu- 
tions of In analogy with the transverse wave case, 
solutions with r = Vq constant exist if x = utt (the vortex 
distance in the x direction is a multiple of the wavelength) . 
In such a case, we can define ujq = a_/(ro^), the charac- 
teristic frequency of the dipole, and V = — 1 + cjq^o, the 



X = V + (-l)^+^ao sinX 



Y = -ujoxot . (25) 



This equation naturally leads to two different cases, de- 
pending on the value of | F/ao | , which determines the ex- 
istence of a fixed point for X. In one case {\V\ > ao\) 
the x-component of the dipole speed X has a constant 
sign, otherwise it changes sign (|V| < ao) and thus can 
eventually approach a fixed point. 

We first notice, that in the very simple case where V = 
(the x-component of speed of the unperturbed dipole is 
equal to the phase velocity), which we call the resonant 
case, the equation (|25]) becomes X = (— l)^+^ao sinX. 
This equation predicts the existence of selected directions 
X = UTT for the dipole propagation, a feature reminiscent 
to the mode locking phenomenon as for instance the one 
observed by in PO) !^ . These directions correspond to the 
case where the dipole is located at an extremum of the 
wave, and therefore where the profile of the wave is flat. 
Therefore, if the condition V = is satisfied the dipole will 
only see this flat profile and the influence of the wave on 
the dipole motion vanishes. Indeed, the general solution 
of dm for V = is 



X = 2 arctan ( e 



tan 



(26) 



and the dipole sets itself in one of the selected directions 
where it does not see the wave. 

A similar qualitative behavior to the special case V = 
0, is found when |F/ao| < 1. The solutions of (|25]) show 
that X(t) converge towards a fixed value 

Xf = (-l)"" arcsin ( — | + mir . 

The difference here is that the dipole is trapped, and the 
wave forces the x-component of the dipole speed to the 
phase velocity uj/k^ by adjusting the slope of the wave 
that the dipole will "see". Fig. [3fa). On the other hand, 
if |y/ao| > 1, the dipole has enough strength to "sur- 
pass" the wave with V as the averaged x-component of its 
speed. Fig. ^h). The wave has only a small influence on 
the dipole. 

In the monopole case it is convenient to introduce polar 
coordinates: x = rcos6>, y = rsin6>, to obtain the system. 



rO 



r = —ao cos X sin(r cos 0) cos , 



ao cos X sin(r cos 0) sin , 



X = — 1 — ao sin X cos(r cos 0) , 
Y = 0. 



(27a) 
(27b) 

(27c) 
(27d) 



These three autonomous nonlinear equations allow chaotic 
behavior of the vortex trajectories, the external wave de- 
stroying the vortex integrals of motion. The variety of pos- 
sible trajectories is illustrated in the series of Figs. HI [5l 
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P^ 0.0 




Fig. 5. (Color online) Vortex trajectory leading to fusion. Typ- 
ical motion in the ao > 1 case; (a) trajectory t G (0,1.5), 
(b) vortex separation as a function of time. Parameters: cci = 
(0,0), X2 = (3,0), r = 87T, To = 4. 



and[6l We introduce various simplifying assumptions in or- 
der to qualitatively describe the behavior of the monopole 
in these different physical regimes. We first investigate the 
long wavelength limit in both weak and finite wave ampli- 
tude approximations, and then we discuss the small wave 
length behavior. We particularly focus on the merging or 
splitting trajectories, in the presence of a finite amplitude 
wave. 

In the weak long- wave regime, we can assume that 
|ao| ~ r <C 1, and |q^| ~ 1; the consistency of these 
assumption will be valid as long as r remain within the 
neighborhood of tq. A series perturbation development 
gives then. 



r = ro , 



X = t — ao cos t , y = Yo , 



(28) 



(29) 



where ujq = o^/ro^ is the natural frequency of the un- 
perturbed monopole. The vortex trajectories are there- 
fore confined in a small region of maximum thickness 2ao 
around the circle. 

In the case where r <C 1 and |ao| ~ \ol\-\ ^ 1 (the wave 
amplitude is finite, comparable to the vortex intensity). 




ln(0 



(b) 



Fig. 6. (Color online) Trajectory of the monopole leading 
to the splitting, (a) Typical splitting process when the wave 
length is smaller than the vortices separation, (b) Long-time 
behavior, the straight line corresponds to a power-law fit with 
exponent 1/3. Parameters are similar to the ones of Fig. (5] but 
with a wave phase of tt. 



the motion equations become. 



(30a) 



2 



-aor cos" cos X ^ — ^^rcosX, (30b) 



X = — 1 — ao sinX , 



(30c) 



where in the distance equation ()30b|) . we approximated 
the cosine term by its mean value, using the fact that in 
the limit r ^ 0, ^ given by (|30ap varies rapidly. We then 
obtain for X an equation similar to ([25|) with V = 1; it 
can be readily integrated to give 



X{t) = —2 arctan ao -\- s tan 
for ao < 1, and for ao > 1, 

X{t) = —2 arctan ao — stanh 



Hit -to))] 



Hit -to)) 



(31) 



(32) 



where s = \aQ — 1|^/^. 

In the case, ao < 1 and r <C 1, quasi-periodic motion 
is possible. Indeed, for weak enough wave amplitude, one 
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retrieves the perturbation result, with X{t) 
t{l — ao)/(l + ao) — 2 arctan(ao) linear and 



ro exp 



0^0(1 + al) 
" 2(1 -ag) 



sinXo(t) 



Xoit) = 



(33) 



corresponding to a periodic separation superposed to a 
constant drift (X(t),0) in the x-direction, as illustrated 
in Fig. [3^. Increasing the initial distance between the vor- 
tices, one can enter a region where the trajectories become 
chaotic (see Fig. [4)3). 

The other case, ao > 1, corresponds to the merging 
of the two vortices, as can be seen in Fig. [5l The mo- 
tion of the center of mass X converges to a fixed value 
Xf = — 2arctan(ao — y^o^—l); as a result, the distance 
r satisfies the long-time solution. 



r ^ To exp 



-y cos(X/)t 



(34) 



where cosXf is positive for ao > 1. Therefore, we obtain 
in this case an exponential approach of the two vortices, 
accompanied with a diverging angular velocity ^ ^ 00 
(Fig.fSb). The merging characteristic time is 2/ao cosX/(ao), 
and within the present approximation, it diverges when 
ao ^ 1. This is consistent with the vortex merging ob- 
served in the numerical solutions, as seen in Fig. [5l 

Let us now consider the situation where, r ^ 1 (the 
wave length is small compared to the size of the monopole), 
^ ^ (^0 ^ OL^ <^ such that the first term in (|27bp may 
be neglected. This case corresponds to the vortex split- 
ting process, shown in Fig.[6l Equations (|27|) become now, 
6 ^ const and 



X = — ao cos X sin X 
X = — 1 — ao sin X cos x 

y 



a+x 



x^ 



r 



(35a) 
(35b) 
(35c) 

(35d) 




similar to Eqs. for the dipole, but with a difference 
in the equation for y{t). Asymptotically (|35a|35bp give 
x{t) -\- X{t) = const., which justifies the approximation 



const , . 

y ^ — 2~' y(^) 



t 



1/3 



(36) 



Indeed, if r is initially large, Eq. ([36|) shows that it remains 
large at long times. We confirmed the validity and verified 
by direct numerical integration of the vortex equations 
([23|) . that the asymptotic behavior of the vortex distance 
follows a power-law with the exponent 1/3 (see Fig. [6)3). 



5 Wave-vortex interaction in the 
Charney-Hasegawa-Mima dynamics 

We consider now the Charney-Hasegawa-Mima equation 
in its two forms, ([9]) and ([TT]) . to study the evolution of 



0.0 5.0 10,0 15.0 20,0 



Fig. 7. (Color online) Potential vorticity of two gaussian 
monopoles at times (left), 7.5 (center) and 15 (right), in 
units of ps and Vd- Without the interaction with a wave, the 
vortices rotate and drift; the small diffusion of the vorticity is 
due to the presence of a viscosity term in Q . 



two vortices in the presence of a transverse and a longitu- 
dinal wave, respectively. We demonstrated, in the context 
of the point vortices model, that the presence of a wave 
can induce a variety of behaviors, including merging, split- 
ting and chaos. However, in real flows or plasma systems, 
a vortex has a spatial extension and waves have to obey 
a specific dispersion relation ([3]). Moreover, as mentioned 
before, the separation of waves and vortices is somewhat 
arbitrary, their interaction modifies both in such a way 
that the distinction can be made only within some length 
scale and time interval. It would then be interesting to an- 
alyze how the presence of a wave might influence processes 
such vortex merging, in the more general framework of the 
Charney-Hasegawa-Mima equation. 
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Fig. 8. (Color online) Potential vorticity of two vortices in the initial field of a transverse wave. The sequence shows 
the monopole complex trajectory, in the chaotic regime of the corresponding point vortex model. From left to right, 
t = 0, 0.8, 1.6, 2.4, 3.2, 4.0, 4.8, 5.6, 6.4, square of size 20^, and vortex size parameter b = 16. 



0.0 5.0 10.0 15.0 20. 



I 



,0 (a) ' 0.0 5.0 10.0 15.0 20.0 ZO.O 0.0 5.0 10.0 15.0 (q) 



0.0 5.0 10.0 15.0 



0.0 (d) 



0.0 5.0 10.0 15.0 



00 (e) 



10.0 15.0 



00 (f) 



Fig. 9. (Color online) Evolution of two gaussian monopoles in the transverse wave case. Vortex and wave parameters are as in 
Fig,[T] The initial conditions, in (a) and (d), differ only in the phase of the wave with respect to the vortex positions. The first 
row (a-c) shows the splitting; the second row (d-f) the merging of the vortices. 



For this purpose we integrate numerically equations 
and (dU), and we kept the dissipation term to ensure 
numerical stability. We consider two identical vortices; for 
each vortex we assign at the initial time, a local vorticity 
f2i in the form of a Gaussian blob, 

7 7-1 

[2i = -V^^ir, 0) + ^{r, 0)/pl = — exp{-br^) , (37) 

TT 

(?/^(r,0) is the initial stream function), where b controls 
the extension of the vortex and F its circulation. In the 
case of a system governed by (|9]), since transverse waves 
are naturally generated, we shall not impose an external 
one, but instead add a wave form to the initial vorticity 
distribution and track its influence on the monopole tra- 
jectories. In the case of longitudinal waves, we solve ([TT]) 
with a source term of the form, 

S{x, t) = So cos(% - ujt) , (38) 



where So measures the source strength (it has dimensions 
of a square frequency). 

Simulations are made using periodic boundary con- 
ditions on a square domain of length L = 20, and the 
number of spectral modes is typically 1024^. The numer- 
ical code is based on a pseudo-spectral scheme with time 
stepping a fourth order Runge-Kutta algorithm; the time 
step is 6t — 10~^, and we used a the numerical viscosity 
of = 0.005. We measure lengths and velocities in units 
of ps and respectively. The initial state consists on 
two equal circulation vortices with vorticity distribution 
([37|) . where 6 = 8 and F = Stt. The set of parameters and 
initial vortex positions, used in the various simulations, 
correspond to the different cases studied with the point 
vortex model. 

The reference case is displayed in Fig. [71 where we 
represent the potential vorticity Oi of Eq. for three 
different times, of two monopoles in the absence of initial 
wave or external source. Without extra perturbations, the 
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two monopoles simply rotate around each other and slowly 
drift in the ^-direction; the distance between the vortices 
remains sensibly constant, their shape evolving in the vis- 
cous time scale. We note that the observed drift motion is 
in fact influenced by the emission of a wake, showing that, 
although Gaussian monopoles are long lived structures, 
they are not exact solutions of the Charney-Hasegawa- 
Mima equation. 



5.1 Transverse wave 

Let us first investigate the transverse wave case, where 
in addition to the two monopoles we superpose an ini- 
tial wave k'^ To cos{ky) . As it could be inferred from the 
large chaotic regions observed in the point vortex model 
(Fig. [1]), for a range of initial conditions the behavior of 
the two monopoles is unpredictable. In the sequence of 
Fig. [HI the vortices rotate, approach each other, inter- 
change vorticity, and separate in the background of the 
wave vorticity. In this case, the initial wave is in phase 
with the vortices, and as a result, they drift with the wave 
and rotate at a high frequency compared to the free case 
of Fig. [71 The horizontal drift is easily explained by the 
model ([T6|) that implies. 



d 

di 
d_ 

di 



X = Foksmuot cos ky^ 



Y = 



(39) 



for the center of mass motion (we use the notation of ([T7|) ). 

The presence of a transverse wave can also trigger pro- 
cesses such as splitting and merging of vortices, as in the 
point vortex model. In Fig. [9] we present numerical simu- 
lations differing in the phase of the initial wave, showing 
the sensitivity on the initial condition typical of a chaotic 
system, and leading to split or fusion. We observe that 
the vortex evolution is much faster than the wave one, as 
can be asserted by comparing the fusion time of tf ~ 3 
in Fig. [9f, to the slow rotation of the two free vortices in 
Fig. [3 

It is interesting to compare more quantitatively the 
motion of the point vortices with the actual evolution 
governed by the full equations. One obvious difference be- 
tween the two approaches is that point vortices lack ex- 
tension, in contrast to the Gaussian vortices of the full 
simulation. However, one may introduce a patch of point 
vortices over a region whose size is of the order of the 
Gaussian monopole, and follow their Hamiltonian dynam- 
ics. Such a comparison is presented in Fig. [TOl where we 
chose the same initial positions and wave phase, and ad- 
justed the point vortices strength to obtain similar time 
scales as in the extended vortex system. We observe that 
the approach of the two coherent structures follows the 
same pattern, including the enhancement of the asymme- 
try between the bottom and top vortices; both patches are 
strongly stretched during they approach, at t ^ 1.5. In- 
deed, in the Charney-Hasegawa-Mima simulation we also 
observe that the top monopole is slightly reinforced as it 



(a) 




(b) 

Fig. 10. (Color online) Comparison of the dynamics of a 
patch of point vortices (a) with the evolution of two Gaus- 
sian monopoles (b), for the merging case at times t = 
0, 0.38, 0.76, 1.13, 1.51. (a) The different patches correspond 
to successive time steps, showing how they approach each other 
from their initial positions (8,10.5), (12,9.5). (b) Successive 
monopole positions. 



approaches the bottom vortex, which moves faster than 
the top one. However a difference arises in the ^-direction 
displacement; while in the vortex model its mean value is 
zero, for extended vortices there is a supplementary drift 
coming from the background vorticity gradient [42 l [43 l [44] . 
Although this effect is responsible of a quantitative dif- 
ference between the point and extended vortex cases, the 
determination of the appropriated range of parameters for 
vortex merging as determined by the point vortex model, 
appears to be rather robust against this effect. 



5.2 Longitudinal wave 

We turn out now to the study of the longitudinal wave 
driven system. The direction of the longitudinal wave can 
be in principle arbitrary. However, it is convenient to well 
separate the longitudinal and transverse situations in or- 
der to identify the physical mechanisms. Therefore, we 
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0.0 5.0 



0.0 5.0 10.0 15.0 



10.0 15.0 



0.0 5.0 10.0 15.0 



00 (c) 



20.0 (f) 

Fig. 11. (Color online) Evolution of two gaussian monopoles with a source of longitudinal waves. Splitting (a-d) or merging 
(d-f) depend on the phase of the source term. (a,d), initial state; (b,e), intermediate state {t = 1.5); (c,f), final state (t = 3). 
Vortex and wave parameters are as in the transverse case (Fig. [9|) . 



take a source that generates longitudinal waves perpen- 
dicular to the propagating intrinsic transverse waves. This 
is the best suited situation to apply the results of the 
point vortex model. Although the source term added to 
the vorticity equation is restricted to a wave of the form 
xsm{kx — uot)^ the evolution of the two equal circulation 
vortices can display a large variety of behaviors. We focus 
here, on the strong interaction of vortices induced by the 
wave, and leading to their separation or fusion as predicted 
by the point vortex model. Therefore, we take parameters 
similar to the ones of Figs. [5] and [6l for the vortices and 
wave intensities, positions and phase. Remarkably, we ob- 
tained the same kind of trajectories in the extended case, 
shown in Fig. [TTl where the vortex pair split or merge 
according to the phase of the longitudinal wave. 

We remark that the relative position of the monopoles 
with respect to the wave influences the behavior of their 
vorticity: the wave possesses its own vorticity which su- 
perposes to the monopoles one, reinforcing the monopole 
in the positive vorticity background and weakening the 
one in the negative vorticity background and causing ad- 
ditional drift [45l[42l[43l[44) . This effect, in the observed 
regime of parameters, slightly affects the motion of the 
vortices during the time of their evolution, without chang- 
ing qualitatively the fusion or splitting process. In ad- 
dition, the vortex interaction mechanisms resulted to be 
much faster than the evolution of the wave vorticity, show- 
ing that neglecting as a first approximation the influence 



of the vortex motion on the wave vorticity appears to be 
justified. 

Although we chose to present results on the monopole 
dynamics in the general case, we note that the point vortex 
model predicts, in the presence of a longitudinal or trans- 
verse wave, that a dipole will follow an almost straight 
trajectory with a superposed oscillation, along certain di- 
rections with respect to the wave propagation. This be- 
havior was also observed in dipole systems obeying the 
Hasegawa-Mima equation . 

6 Conclusions 

In this paper we introduced a simple model based on the 
dynamics of point vortices, in order to get an idea on 
how extended vortices and waves interact in more real- 
istic flows. As a matter of fact, it turned out that this 
simple model was much richer in the variety of behaviors 
it could reflect, than what would be expected. Namely, 
it demonstrated how stable a dipole was, and in oppo- 
sition how like-circulation vortices could be destroyed or 
merged by a wave. Comparison of the point vortex plus 
wave analytical results and the Charney-Hasegawa-Mima 
simulations, shown that the distinction between localized 
and extended vorticity is justified by the existence of well 
separated length and time scales, and is also relevant to 
determine, for instance, the conditions for vortex merging 
or splitting. 
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We considered in particular, the interaction of equal 
circulation vortices in the presence of an intrinsic drift 
(transverse) and of an external (longitudinal) wave. We 
found that, depending on the relative phase of the wave 
with respect to the positions of the monopoles, merg- 
ing, splitting or even chaotic trajectories are possible. Us- 
ing the same set of parameters, we confirmed that the 
point vortex model reasonably predicts the behavior of the 
extended vortices obeying the Charney-Hasegawa-Mima 
equation. 

In summary, the results obtained in the case of the 
point vortex model and confirmed by numerical simula- 
tions, could partially explain how coherent structures can 
appear by the merging of small structures and then or- 
ganize themselves into dipoles and monopoles, especially 
in systems where waves can be present or generated by 
the dynamics of the structures. We may hope that this 
simple model, that permits to identify the range of pa- 
rameters relevant for the vortex- wave interactions, could 
be a help in understanding the dynamics of large coherent 
structures in decaying turbulent flows. 
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where po = const. The equation of state, for instance the 
adiabatic relation between pressure and density P ^ p^ ^ 
complete the basic equations. Applying the rotational to 
(HO]), one obtains the vorticity evolution equation. 



d 
dt 



VP 

P 



gz^ fv X z, 



(40) 



where 



d_ 

di 



d_ 

dt 



is the total derivative, p is the fluid density, P the pressure, 
g the acceleration of gravity, and / = f{y) the Coriolis 
frequency (which is a function of the latitude). The usual 
choice of coordinates is x in the equator east direction and 
y in the north direction. Note that the plasma convention, 
as in (j2j), exchange the role of the coordinates: (x^y) 
(?/, —x). The continuity equation writes. 



V • V 



dt Po 



(41) 



d_ 

dt 



{lv + fz) + (cj + fz)V • V- 
(cj + /£) • Vv = - 



-V-xVP, (42) 
P 



where = V x v. 

We consider a fluid layer at height Hq of width h = 
h{x^t)^ with X = ix^y). To describe the dynamics of this 
layer it is convenient to separate in the equation of motion 
(|^Q|) . the horizontal (±) and vertical [z) components: 



d 1_ ^ . . 

—v±_ — V±P + X z, 
dt p 



and 



dvz 
~dt 



IdP 

p dz 



(43) 



(44) 



respectively. We also consider that in addition to the hy- 
drostatic pressure p^gh^ there can be a forcing pressure 
p = p{x, t) such that the total pressure writes. 



P = Po ^ pogh ^ p . 
The extra pressure term induces a vertical flow, 
dvz 1 dp 



dt 



Po dz ' 



(45) 



(46) 



7 Appendix: Charney equation with source 
term 

In this appendix we briefly discuss the derivation of the 
Charney-Hasegawa-Mima equation in the geostrophic frame- 
work [21], in order to establish its relation with (|8]) and 
the origin of the source term. 

The velocity field v satisfies the motion equation. 



where we neglected the nonlinear convective term and re- 
placed everywhere the density by its mean value po . More- 
over, from the horizontal motion equation (|^3|) . the bal- 
ance of the hydrostatic pressure gradient and the Coriolis 
force leads to a drift velocity. 



g_ 

fo 



z X \/h , 



(47) 



where, to this order, the Coriolis frequency is replaced 
by its averaged value /o, and the inertia and unsteady 
pressure p terms, are neglected [53]. The total velocity 
field can thus be written as, 



V = v± 



(48) 



where the second term, that takes into account the fluctu- 
ations around the equilibrium drift velocity (|T7|) . is of the 
same order as the fluctuating pressure term in (|^5]) . Equa- 
tion (|48|) can be compared with (|8|), to note that the point 
vortex model reduces to the advection of the localized part 
of the height h field, in the field of the drift velocity 
possibly driven in addition, by a external velocity field 
V^. Using this expression in (|36]), one obtains. 



d 

di' 



po 



(49) 



which relates the pressure fluctuation with a vertical (po- 
tential) motion. 
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Introducing (|Trj) into (|^ . and neglecting the baro- 
clinic term (which vanishes for adiabatic motions), the 
vorticity equation becomes, 



dt 



(50) 



where uz = \/ x v±. Using now (jH]), with 



ln^«ilnfl + A + iL 

^0 



Po 7 



neglecting uj in the second term of (|5Q]) , and averaging 
over the vertical layer, one obtains. 



d_ 

dt 



f 



h 



dvz 
dz 



(51) 



where all the quantities depend on the horizontal coordi- 
nates. The right-hand side term can be explicitly calcu- 
lated using (gni, 

d p dvA_-l/d'^4> 2d'^<t>\ _ S{x,t) 
dt^po^ dz/~ cl\df^ "^9^2/- /o ' 

(52) 

where Cs = {jPo/ po)^^'^ is the sound velocity, and where 
we defined S = S{x^t)^ the effective source of vorticity 
that contains the contribution of vertical unsteady mo- 
tions. The appearance of the sound speed is a consequence 
of the assumption on the equation of state. An equivalent 
equation, with a different source term, can be obtained if 
instead of an external force, one introduces in the energy 
equation a heat source, and relates the pressure fluctua- 
tions with the temperature fluctuations PT | [M ) . It is worth 
noting that if (j) satisfy the wave equation, the source term 
reduces to the form 5* ^ k'^cj), where k = (kx^ky) only 
contains the perpendicular component of the wave vector. 
Although the specific form of the source term will depend 
on the physical mechanism from which it originates, being 
proportional to dvz/dz^ it naturally generate longitudinal 
perturbations. 

Replacing ou = (^//o)V^/i and (|52]) into ([5T]) . and 
putting / = /o + we finally get 



dt V/o 



7^0 ' 



h 



S{x,t). (53) 



In the usual quasi-geostrophic approximation, the S term 
is neglected, and one recovers the Charney-Hasegawa-Mima 
equation ([9]), where we used the following identifications 
between plasma and fluid parameters: i/j = gh/ fo, Vd = 
PjHq/ fo, and = jHq/ fo^ and changed axes according 
to the plasma convention, (x^y) {y, —x). 
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